BCS Superconductivity in Heavily Boron-Doped Silicon 
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Scanning tunneling spectroscopies in the subKelvin temperature range were performed on super- 
conducting Silicon epilayers doped with Boron in the atomic percent range. The resulting local 
difTerential conductance behaved as expected for a homogeneous superconductor, with an energy 
gap dispersion below ±10%. The spectral shape, the amplitude and temperature dependence of the 
superconductivity gap follow the BCS model, bringing further support to the hypothesis of a hole 
pairing mechanism mediated by phonons in the weak coupling limit. 
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Upon heavy Boron doping, Silicon shows a supercon- 
ductive resistive and diamagnctic transition at low tem- 
perature [1]. This recent discovery has been held as a 
significant milestone [2] because of the archetypical na- 
ture of this covalent cubic semiconductor [3-6]. Simi- 
larly to diamond, the charge carriers involved in Silicon 
superconductivity are generally thought to be holes at 
the Fermi level lying within and near the top of the va- 
lence band [7, 8]. Unlike diamond [9], the superconduc- 
tive transition in Silicon requires a Boron content (sev- 
eral percents [10]) much higher than the one needed for 
the insulator to metal transition (about 80 ppm [11]). 
This threshold also lies above the Boron solubility limit 
in crystalline Silicon, raising the question of the homo- 
geneity of doping and hence superconductivity in Si:B 
epilayers obtained by out-of-equilibrium growth. On a 
more fundamental basis, the nature of the pairing mech- 
anism in this new superconductor remains to be deter- 
mined. So far, the critical temperature values, the opti- 
cal phonon softening observed by Raman spectroscopy in 
the normal state [1] and the ab initio calculations of the 
electron-phonon coupling [7] have only shown that the 
overall features of this new superconductor were com- 
patible with a conventional pairing mechanism. These 
open questions call for more local investigations of su- 
perconductivity in doped Silicon. 

In this Letter, we report a detailed local study in heav- 
ily Boron-doped Silicon by subKelvin tunneling spec- 
troscopy. The local electronic density of states (LDOS) 
spectral shape, the amplitude and temperature depen- 
dence of the related energy gap in the excitation spec- 
trum are consistent with the picture of a superconductor 
following the BCS model. The superconducting energy 
gap was found to vary locally by up to ±10%. 

Superconductive Silicon samples were prepared by Gas 
Immersion Laser Doping (GILD) [12]. A precursor gas 
atmosphere BCI3 was injected on a {100}-oriented Silicon 
wafer surface in order to saturate Silicon chemisorption 
sites by Boron atoms. After each injection, the surface 



was melted by an ultraviolet laser pulse, which inserted 
Boron atoms by diffusion into the Silicon wafer. This out- 
of-equilibrium method allows substitutional Boron atoms 
to be incorporated into the Silicon lattice well beyond the 
thermodynamic solubility limit of about a few 10^° cm~^ 
(i.e. below one atomic percent). The doped area deter- 
mined by the laser spot size was typically 2 by 2 square 
millimeters. As a result, high concentration Secondary 
Ion Mass Spectroscopy (SIMS) profiles [10, 13] show, be- 
low a few nanometer-thick surface oxide, an homogeneous 
doped layer thickness and a concentration plateau over 
30 to 130 nanometers. 

In the present work, we studied different samples (Al, 
B{1,2,3} and CI), which were prepared by respectively 
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FIG. 1. (Color online) Bottom: normalized differential con- 
ductances {dl /dV) of tunnel contacts on superconducting Sil- 
icon (sample Bl, Tgampie = HO mK) at one micron apart po- 
sitions: a and b. The dotted line is a BCS fit. Top: variation 
of the dl /dV spectrum with respect to sample temperature 
from bottom to top Trample = HO, 350, 425, 440, 460, 500 mK 
on the same sample and at position a. The set-point tunnel 
resistance was 1 MH. 



200, 500 and 400 laser shots, yielding different 'macro- 
scopic' critical temperatures Tc— 350, 510 and 590 mK. 
The above-mentioned values of the critical temperature 
Tc were determined as the temperature corresponding to 
a resistance below 10 % of the normal state value. This 
criterion was chosen because this critical temperature co- 
incides with that of the magnetic susceptibility onset [10]. 

In sample Al, which corresponds to the first set of su- 
perconducting Silicon studied in Ref. [1], the laser melt- 
ing duration was short (25 ns) and the layer very thin 
(30 nm) , giving a low critical temperature T^ value and a 
large transition width of AT « 80 mK. This broadening 
is defined from 10% to 90% of the resistive transition. 
By increasing the number of laser shots, it was possi- 
ble to increase the layer thickness (80 nm) and optimize 
the critical temperature [10]. For instance, sample CI 
has the highest critical temperature Tc= 590 mK and 
the sharpest transition: AT w 20 mK. In this sample, 
a mean free path I of 3 nm and a superconducting co- 
herence length ^5 of about 50 nm have been previously 
estimated via the square resistance i?n ~ 17 J7 at room 
temperature and the upper critical magnetic field, re- 
spectively [10]. In the case of a large number of laser 
shots, a fraction of Boron atoms is not incorporated on 
substitutional sites [10]. For samples B, which exhibit a 
130 nm thickness, this leads to a more disordered Sili- 
con structure, a lower critical temperature Tc value and 
a larger transition width: AT « 40 mK. 

From the critical temperature values discussed above, 
the energy of the local superconducting gap is expected 
to be in the low milli-electronvolt range. The use of a 
scanning tunneling microscope (STM) working in a dilu- 
tion refrigerator [14] is therefore crucial to achieve a lo- 
cal spectroscopy with the proper energy resolution. This 
surface probe has been for instance applied to supercon- 
ducting diamond [15, 16] and it is particularly well suited 
to study thin films such as Si:B epilayers. Furthermore, 
STM probing has the advantage to give access to the 
intrinsic sample inhomogeneity. 

For local tunneling spectroscopy, a sweep of bias volt- 
age V was applied between the Silicon sample and a 
metallic tip [17] at a fixed sample-tip distance, and the in- 
duced tunnel current I was recorded. We used two types 
of measurement that gave similar results: an a.c. lock-in 
technique with a 15 /iV excitation at about 2 kHz fre- 
quency and a numerical derivation from a d.c. I{V) mea- 
surement. At thermal equilibrium, the local differential 
conductance (dl/dV) gives access to the local density 
of states p smeared out by the thermal energy 3.5 keT: 
dl/dV oc /jjdE9f/9E • /9(E - eV). In this expression, 
the energy distribution function derivative with respect 
to the energy di/dE features the thermal broadening. 
For very low temperature experiments, we define an ef- 
fective temperature in order to include several additional 
effects, which also smear out the electronic current char- 
acteristic. These contributions are related to local heat 



induced by out of equilibrium electron injection, elec- 
tromagnetic noise (black-body radiation of hot objects, 
imperfect filtering) , weak electron- phonon coupling. The 
effective temperature is usually higher than the refrig- 
erator base temperature and it determines the practical 
energy resolution of the measurement. In our case, this 
resolution is better than 2kBToff(300 mK) « 50 ^eV. 

Fig. 1 bottom presents the local differential 
conductance measured at a sample temperature 
Tsampic = 110 mK on sample Bl (i.e. sample numbered 1 
in the series B). A superconducting energy gap is clearly 
visible at low bias voltage with two coherence peaks 
around \V\ = 0.13 mV. A constant differential conduc- 
tance characteristic of an ohmic behavior is recovered 
at higher voltage. Such a spectroscopy measurement 
represents the first local evidence of superconductivity 
in heavily Boron-doped Silicon. A series of positions 
were probed. Fig. 1 displays as an illustrative example 
spectra acquired at two positions a and b distant from 
one micrometer, which is one order of magnitude larger 
than the coherence length ^5 . Similar spectra were found 
in all cases, which indicates that the superconductivity 
is spread over the whole sample with a similar behavior. 
The signature of inelastic tunneling processes involving 
softened optical phonons at 450 and 600 cm~^ [1] (i.e. 
55-75 meV) has also been sought but without success. 

The differential conductance dl/dV spectrum can be 
fitted by a BCS-type density of states (dotted line in 
Fig. 1). This expression does not include any Dynes 
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FIG. 2. (Color online) Left axis: variation of the supercon- 
ducting energy gap versus the sample temperature, recorded 
for three different positions in sample B3. Position 1 and 2 
are 15 nm apart, whereas position 3 is 100 nm away. The 
continuous line is to a fit to the BCS theory. Right axis: elec- 
trical resistance of the same sample versus temperature (dot- 
ted line), showing a macroscopic superconductive transition 
with a critical temperature (see text for its precise definition) 
Tc ^ 510 mK. 



parameter but takes into account a fixed effective elec- 
tronic temperature of 300 mK. A gap amplitude about 
AlT'sampic = 110 mK) = 87±7 fieV is extracted from the 
fit. The error bar in the gap value is significantly smaller 
than could be naively expected from the thermal broad- 
ening, since the fit is obtained with the gap as the single 
adjustable parameter. 

Recent tunnel junction experiments performed on spe- 
cific arrays of microdefects in doped Silicon have yielded 
a large energy gap, which was attributed to a high- 
temperature superconductivity [18]. Here, in order to 
confirm the superconductivity origin of the measured 
gap of excitations, we varied the sample temperature 
and recorded the corresponding local differential conduc- 
tance. A typical set of differential conductance data at 
different temperatures, obtained for sample Bl at a single 
position a, is presented at the top of Fig. 1. As expected, 
the amplitude of the energy gap decreases by increasing 
the sample temperature (lower to upper curve at zero 
bias), while the smearing also increases. 

At every sample temperature, the value of supercon- 
ductivity gap can be obtained by fitting the spectrum 
with a BCS expression. This allowed us to determine 
the temperature dependence of the superconducting gap. 
Fig. 2 presents this data for three different positions (1, 
2 and 3) in sample B3. The distance between 1 and 2 on 
one hand and (1,2) and 3 on the other hand was respec- 
tively below and above the estimated superconducting 
coherence length. At these three different positions, the 
energy gap values coincide within the error bars. 

The temperature dependence of the gap is also well fit- 
ted by the thermally-smeared BCS prediction displayed 
as the continuous line in Fig. 2. It follows for sample B3 
a zero temperature energy gap A(0) = 78 ± 9 ^eV and a 
local critical temperature Tc = 480 ± 30 mK. The criti- 
cal temperature accessed in this way is a local quantity. 
It corresponds to a value close to the tail of the super- 
conductive resistive bulk transition (see dotted curve in 
Fig. 2), which sits at Tc ~ 510 mK. This distinction 
between the 'resistive' and the local critical tempera- 
ture properly refiects the difference between the trans- 
port measurement and the local probing. 

In total, we have measured five samples, and more 
than one hundred spectra have been obtained on dif- 
ferent spatial positions. All curves showed similar re- 
sults, demonstrating that the data displayed in the previ- 
ous Figures are characteristic of superconductive heavily 
Boron-doped Silicon. The quality of the tunnel junction 
has been controlled by recording the differential conduc- 
tance dl/dV at different tunneling resistances, i.e. differ- 
ent sample-tip distances. When changing the set-point 
tunnel resistance between 10 MJ7 and 1 Mfi, the tun- 
neling spectra overlap after renormalization with respect 
to the normal differential conductance, as expected. To 
warrant a clean surface, buffered oxide etch (BOE) or 
diluted hydrofiuoric acid (HF) cleaning has sometimes 



been performed but no significant difference has been 
revealed compared to the 'as grown' surface. Further- 
more, we succeeded in fabricating a millimeter square 
solid superconductive-normal metal tunnel junction to 
probe the tunneling current averaged on a large scale. 
The measured current-voltage characteristic provided a 
similar differential conductance spectrum with a super- 
conducting bulk gap value comparable to the one mea- 
sured locally by STM [19]. 

Finally, the difference in critical temperature of the 
three series A, B and C enabled us to investigate the evo- 
lution of the zero temperature superconducting energy 
gap value with the critical temperature. Fig. 3 displays 
the series of (A(0), Tc) couples for the whole experimen- 
tal runs set on the five investigated samples. The BCS 
model prediction A(0)/fc_B Tc — 1.76 is indicated as a full 
line. The critical temperatures Tc were determined either 
by macroscopic transport measurement (empty symbols 
in Fig. 3) or, when available, by local temperature de- 
pendence measurements, following a similar procedure 
such as the one shown in Fig. 2 (full symbols in Fig. 3). 
All measured data points are located around the full line 
corresponding to the BCS model. Data with a critical 
temperature measured locally agree with the BCS ratio 
within the error bars. For the rest of the data. Gap values 
smaller than the BCS prediction are ascribed to a weak- 
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FIG. 3. (Color online) Dependence of the extrapolated zero- 
temperature superconducting gap A(0) with the critical tem- 
perature Tc for series A, B and C samples. For full symbols, 
Tc was locally determined through the local energy gap tem- 
perature dependence. For empty symbols, Tc was obtained 
via macroscopic transport measurements. The straight line 
indicates the BCS ratio A{0)/kBTc = 1.76, while the dotted 
lines show a ± 20% deviation. Sample Bl and CI data with 
empty (O, V) and full (•) symbols are measured at a macro- 
scopic distance (i.e. at a different cryostat cooling). Sample 
B2 data (A) were measured 0.5 /^m apart. Sample B3 sym- 
bols (■) summarize Fig. 2 results. 



ening of superconducting energy gap, which may occur 
locally because of a higher disorder or a lower charge 
doping at the surface compared to the bulk of the film. 

On a given sample, the measured energy gaps show 
a maximum ±10 % scattering. The same scattering is 
observed over small and large distances, see for instance 
sample Bl data in Fig. 3. The disorder accessed through 
the local superconductive energy gap thus sits at a very 
small scale, of the order of the mean free path. The en- 
ergy gap variation normalized to the mean value is larger 
than the width of the resistive transition normalized to 
the critical temperature. Actually, the local energy gap 
dispersion measured by STM gives an accurate picture 
of the actual spatial variation of the superconductivity, 
whereas the resistive transition indicates only the appear- 
ance of a superconducting percolating path throughout 
the sample. These two approaches can give a different 
picture, depending on the investigated sample. For in- 
stance, sample CI yielded the sharpest superconducting 
resistive transition [10] of all samples but an energy gap 
dispersion similar to the one observed in series B sample 
with a broader resistive transition. Finally, let us point 
out that although the gap inhomogeneity observed here is 
reminiscent of predictions on superconductivity in highly 
disordered metals close to the metallic to insulator tran- 
sition [20] , our samples have a doping level two orders of 
magnitude higher than the relevant critical threshold. 

In conclusion, the tunnel differential conductance spec- 
tral shape, the BCS-compatible dependence of the su- 
perconducting gap with the temperature, as well as the 
A(0)/fcB Tc ratio close to the BCS value, are three strong 
experimental indications that superconductivity in heav- 
ily Boron-doped Silicon follows the weak coupling BCS 
model. The microscopic gap variation of a few percents 
fits well with expectations for a strongly disordered metal 
like highly Boron-doped Silicon. This progress in under- 
standing superconductivity in Silicon opens the way to 
the future completion of reliable monolithic devices based 
on hybrid Superconductor/Normal junctions with heav- 
ily Boron-doped regions as the superconductor. 
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